NCU-HEP-k033 
larXiv: 0904 .43521 [hep-ph] 
April 2009 

Electric and anomalous magnetic dipole 
moments of the muon in the MSSM 

Kingman Cheung"'^''^, Otto C. W. Kong'^, and Jae Sik Lee^ 

° Department of Physics, National Tsing Hua University, Hsinchu, Taiwan 300 

'^Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan 300 

'^Division of Quantum Phases & Devices, Konkuk University, Seoul 143-701, Korea 

'^Department of Physics and Center for Mathematics and Theoretical Physics, 
National Central University, Chung-Li, Taiwan 32054 

ABSTRACT 

We study the electric dipole moment (EDM) and the anomalous magnetic dipole moment 
(MDM) of the muon in the CP-violating Minimal Supersymmetric extension of the 
Standard Model (MSSM). We take into account the contributions from the chargino- and 
neutralino-mediated one-loop graphs and the dominant two-loop Higgs-mediated Barr-Zee 
diagrams. We improve earlier calculations by incorporating CP-violating Higgs-boson 
mixing effects and the resummed threshold corrections to the Yukawa couplings of the 
charged leptons as well as that of the bottom quark. The analytic correlation between the 
muon EDM and MDM is explicitly presented at one- and two-loop levels and, through 
several numerical examples, we illustrate its dependence on the source of the dominant 
contributions. We have implemented the analytic expressions for the muon EDM and 
MDM in an updated version of the public code CPsuperH2 . 0. 
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1 Introduction 



The anomalous magnetic dipole moment of the muon, , has provided one of the most 
sensitive test grounds for the vahdity of the Standard Model (SM) [1]: 

af^^ = 11 659 208 (6.3) x 10"^° . (1) 

At the same time, it also provides an important constraint on new physics with precise SM 
predictions available. The SM prediction consists of QED, electroweak (EW), and hadronic 
contributions. The hadronic contribution is further decomposed into leading-order (LO) 
part and higher-order vacuum polarization (VP) and light-by-light (LBL) parts [2-18] B 

^SM _ ^QED , ^EW I ^Had.(LO) , ^Had.(VP) , ^Had.{LBL) 

= 11659177.3(5.3) X 10~^°. (2) 

Equations ([1]) and ([2]) suggest that there is currently a 3.7 cr discrepancy between the ex- 
perimental result and the SMprediction, which can be attributed to possible contributions 
from physics beyond the SM [J: 

Aaj^^P = af"^ - af" = 30.7(8.2) x 10"^° (3.7a) . (3) 

One of the most appealing scenarios for physics beyond the SM is augmented with 
a softly broken supersymmetry (SUSY) around the TeV scale. The supersymmetric con- 
tributions to from such models are known up to dominant two-loop contributions. The 
one-loop results can be found in [19-24] and the two-loop results in [25-30]. It is well- 
known that the dominant two-loop contribution comes from Higgs-mediated Barr-Zee di- 
agrams [31]. The error associated with the known SUSY contributions is estimated to be 
~ 2.5 X 10~^° [32], which is smaller than half of the current experimental and SM the- 
oretical ones. On the other hand, the SUSY augmented models can contain additional 
CP-violating phases beyond the SM Cabibbo-Kobayashi-Maskawa (CKM) phase leading 
to sizable EDMs [33-36]. The current limit on the muon EDM [37] is 

< 1.8 X 10"^^ ecm (95%C.L.), (4) 

which is much weaker than the constraints from the non-observation of the Thallium [38], 
neutron [39], and Mercury [40] EDMs: 



(5) 



Mti| 


< 


9 X 10" 


-25 


e cm 


\dn\ 


< 


3 X 10" 


-26 


e cm 




< 


2 X 10" 


-28 


e cm 



*See Tabic [T] for details of the SM prediction. 

^ The number incorporates the hadronic VP result calculated based on measurements at electron- 
position storage ring. Calculation based on hadronic r-decays gives a different result, which is considered 
less reliable. If the r-based result is used, the overall discrepancy reduces to only about la. 
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Table 1: The SM prediction of a^: see, for example, Ref. [2]. For the quoted value in 
Eq.([2]) we use the results in [8] and [17] for the hadronic leading-order and the hadronic 
light-by-light contributions, respectively. 





X 10^° 


Sttf, X 10^° 


Ref. 


„QED 


11658471.810 


(0.016) 


[3-5] 




15.4 


(0.2) 


[6] 


„Had.(LO) 
Had.(LO) 


690.9 
689.4 


(4.4) 
(4.6) 


[7] 

[8] (*) 

L J V / 


„Had.(LO) 
"/J 


692.1 


(5.6) 


[9] 


„Had.(LO) 


694.4 


(4.9) 


[10] 


^Had.(LO) 


691.04 


(5.29) 


[11] 


„Had.(VP) 
"m 


-9.8 


(0.1) 


[8,12] 


„Had.(LBL) 
^Had.(LBL) 
^Had.(LBL) 
^Had.(LBL) 
„Had.(LBL) 
„Had.(LBL) 

"a* 


8.0 
13.6 
11.0 
11.6 
10.5 
10.2 


(4.0) 
(2.5) 
(4.0) 
(4.0) 
(2.6) 
(~3) 


[13] 
[14] 
[15] 
[16] 

[17] (*) 
[18] 



Nevertheless, if the muon EDM experiment in the future can achieve the projected sensi- 
tivity [41] 

rf^ ~ 1 X 10-2^ ecm, (6) 

the precision of the experiment will be comparable to that of the current Thallium EDM 
experiment. 

In this paper, we study the correlation between the muon EDM and MDM in the 
CP-violating MSSM [30,42]. We present the relation between the one-loop chargino- and 
neutralino-mediated EDM and MDM of the muon. We also derive an analytic relation 
between them in the two-loop contributions from the dominant Higgs-mediated Barr-Zee 
diagrams. We improve the earlier results by including CP-violating Higgs-boson mixing 
effects in the Barr-Zee diagrams and resumming the threshold corrections to the muon, 
tau, and bottom-quark Yukawa couplings in the one- and two-loop graphs. We then focus 
our numerical studies on three types of scenarios in which (i) the muon EDM and MDM 
are dominated by the one-loop contributions, (ii) the lightest Higgs boson is mostly CP 
odd and lighter than ~ 50 GeV, and {in) the dominant contributions to the muon EDM 
and MDM come from the two-loop Barr-Zee graphs. 

For the presentation of our analytic results, we follow the conventions and notations 
of CPsuperH [43] , especially for the masses and mixing matrices of the neutral Higgs bosons 
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and SUSY particles. The layout of the paper is as follows. Section 2 presents formulas 
relevant to the one-loop contributions to the muon EDM and MDM from chargino- and 
neutralino-mediated diagrams. Non-holomorphic threshold effects on the muon Yukawa 
coupling have been appropriately resummed. In Section 3, we present analytic results for 
the Higgs-mediated two-loop Barr-Zee diagrams. For this, most importantly, the resummed 
threshold corrections to the tau and bottom-quark Yukawa couplings and the CP- violating 
Higgs-boson mixing effects have been incorporated. In Section 4, we present some numerical 
examples, depending on the source of dominant contributions to the muon EDM and MDM. 
We summarize our findings in Section 5. 



2 One-Loop EDMs and MDMs of charged leptons 



The relevant interaction Lagrangian of the spin-1/2 lepton with EDM di and MDM ai is 
given by 

(7) 



£spin-i/2 = -\di {la,, 75 I) + 1^ {la,, I) F^^ 



where cx^'' = ^[Y .Y] = i{l''Y - 9^")- The EDM and MDM amplitudes are given by 



{M 



MDM J 



eg; ^ 
2mi 



u{p'){ia^''q,)u{p) , (TUedm)'^ = diu{p'){a^''^,q,)u{p) . 



Generic interactions of charginos Xi2 or neutralinos X? 2 3 4' collectively denoted by with 
a lepton / and a slepton (sneutrino) l'^ 2 {I'l) are given by 



The Lagrangian for the interactions of x and 2 with the photon field A, is 

£xxA = -eQxfex)^'' and Cpp^ = -ieQpl* ^^^l'A'' . 
The diagrams in Fig. [T] induce the MDM and EDM of the lepton / as follows: 



mr 



9Rij 



9Lij 



-Q^A{mlJml) + QpB{mlJm};, 



+ 



m 



Xi 



mi 



3fte 



'diY 



m 



^^3m 



9Rij] 9Lij 



xii' \ xii' 

9Rij j 9Lij 



Qx M^lJrnj) + Qp B{mlJm}, 
Qx MmlJml) + Qp B{mlJmf, ] 



where 



A{r) 



2(1 - r) 



3 -r + 



2 Inr ' 
1 — r 



(9) 



(10) 




Figure 1: Feynman diagrams for the EDM {di}^ and MDM (ai)^ of the lepton I induced by 
X- The photon line can be attached to the slepton I' line or the x line. 

^, , 1 2rlnr\ 
Bir) = — 1+rH , 

^ ' 2{i-ry\ 1-r y ' 

A/ \ 1 ( ^ 2 6rlnr\ 

■^('■' = I^(r37)5(2 + 5r-r^ + — j. 

= Ii(r37)j(l-5r-2r'-— j. (12) 

We have checked that our analytic expressions for the one-loop MDM and EDM agree with 
those given in, for example, Ref. [23] and [36] with Q^- = Qj- = —1 and Q^o^ = Qy^ = 0. 
Note A{1) = -1/3, 5(1) = 1/6 and ^(1) = B{1) = 1/24. Finally, in Eq. ([U]), the 
chargino-lepton-sneutrino couplings are given by 

9lT = -9iCRh, 9lT = hnCLh, (13) 
and the neutralino-lepton-slepton couplings are given by 

gf^ = -V2gTiN:^{U%-V2gtw{Qi-Ti)N:,{U%-hiN^^^^^^^ 

^?gf = V2gtwQiNa{U%-h:N,,{U%, (14) 

with = -1/2 and Qi = -1. 



There are non-holomorphic threshold corrections to the Yukawa couplings hi which 
appear in the chargino and neutralino couplings [44-46]. These corrections become sig- 
nificant at large tan/3 and we resum these effects by redefining the Yukawa couplings as 
follows: 

hi = y^^^. (15) 

In the presence of the CP phases. A; takes the form [47] 



A, 



Of. 



w 



I Mo 



I Mo 
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+ 2 










1 

+ 2 











+ 



J(m^,m^^,|Mi| 



l42l'A^|,|Mir,|/xp) 



where the one-loop function is 



I{x,y,z) = 



xy \n{x/y) + yz ln(y/2;) + xz ln(z/x) 
(x -y){y- z) {x - z) 



(16) 



(17) 



These effects on the muon MDM have been considered in Ref. [29] only in the CP conserving 
case. 



3 Barr-Zee Graphs 

The dominant two-loop contributions come from the Barr-Zee diagrams mediated by neu- 
tral Higgs-boson exchanges, see Fig. [2l Here, we consider loops of third-generation fermions 
and sfermions, charged Higgs bosons, and charginos. The diagrams can be evaluated first 
by one-loop computation of the Hi{k — q)-^'^{q)-'y'^{k) vertex. Thanks to gauge invariance 
the effective vertex takes the form 

rr = [g'^k ■ q - k^q""] S.ik') + [e^^-^k^qp] V^ik') , (18) 

where q is the incoming four-momentum of the external photon, p [p') the four-momentum 
of the incoming (outgoing) lepton, and k the four-momentum of the internal photon going 
out of the upper loop. Note that p' = p + q. Explicitly, keeping only the terms linear in 
the external momenta, the fermionic contributions are given by 

vm ^ -i-X''Q)^,^,9knj i,^,J),._.^y (19) 

with the generic f-f-Hi interaction: Cjjjf = —gf f {gfjjj+ilb 9h ff) f • '^^^ color factor 
Nfj = 3 for quarks and 1 for leptons and charginos. On the other hand, the sfermion loops 
contribute only to the scalar form factor as 

- ^x^^QV^^«.H [ ■ (20) 

with the generic f-f-Hi interaction: jCj^j^j = v gj^j,j f* f Hi . 
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%p) Kp') 






Figure 2: Barr-Zee diagrams: the Hi lines denote all three neutral Higgs bosons, includ- 
ing CP-violating Higgs-boson mixing, and heavy dots indicate resummation of threshold 
corrections to the corresponding Yukawa couplings. 



Together with the interaction Lagrangian for the couphngs of the neutral-Higgs bosons 
and photon to charged leptons 



^Hji = -91 1 {9hJi + «75 9hJi) I Hi, C 



All 



-eQih^lA 



(21) 



where gi 
iM^ 



g nil 
2Mw 

BaiT— Zee 



[as in the CPsuperH convention), the Barr-Zee amphtude is given by 



(22) 



egiQ 



d^k u{p') [rr7. (-^) {9h, + ^75 9'nji) + ^9^1 + 9^^) (^) ^Tl. \ n{p) 



F F (F - M%. 



where again we keep only the terms linear in the external momenta. We note the numerator 
of the integrand is proportional to 



i9h S^ - 9hJi ^0 i^^'^Q'^) - i9Hji S^ + 9%ii V.){a^^imu) . 



(23) 



We observe that the first term gives the MDM while the second one gives the EDM of the 
lepton /. Consequently, the Barr-Zee contributions to the MDM and EDM are related by 
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H 



dUiU 9Hill 

9hJi ~^ ~9Hiii 



(24) 



with normalizations of the MDM and EDM of the lepton / as given by Eq. ([7]). For 
the dipole moment diagrams with chirahty flip inside the loop, the MDM and EDM parts 
correspond directly to real and imaginary parts of the overall amplitude. The above relation 
is a recasting of that statement in terms of the effective scalar and pseudoscalar couplings 
of the generally CP-mixed Higgs states involved. To give the MDM result explicitly, with 
Qi = —1, we have § 
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q=t,b 



+ 
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9Hil+l- 
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EE 
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(25) 



where Txi = m'^/Mjj_ and the two-loop functions F{t), /(t), and g{T) are 
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(26) 



For genuine SUSY contributions, the embedded SM contribution should be sub- 
tracted. We estimate the SM contributions as 



SM 



E 

q=t,b 



em '^q ■■"l r , \ 
fiTqSM) 



fir. 



tSM) 



(27) 



■i-Here we add the contribution from the charged Higgs boson loop and confirm the positive signs of the 
fermionic Barr-Zee contributions [36]. 
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where t/sm = ^}/^Hsm- ^^^^ 



^H 
SM 



1.8 X 10-^1 when Mhsm = 100 GeV and 



it decreases as Mh^m increases. In our numerical analysis, we safely neglect (ai)f^^. 



4 Numerical Analysis 

Adding up all the contributions considered in the previous sections and neglecting {ai 
the supersymmetric contribution to the muon MDM is given by [32] 



1^ 



/X A 4aem 1 MsuSY 

(s) SUSY = 1 :— log 



[a,Y^ + [ay + (a^)^ 



('^m) SUSY + ('^m) SUSY + ('^At)fuSY i (2^ 



where the large QED logarithm takes into account the renormalization-group (RG) evolu- 
tion of from the SUSY scale down to the muon-mass scale. The logarithmic correction 
amounts to -7 % and -9 % for Msusy = 100 GeV and 1000 GeV, respectively. 



4.1 A typical scenario where 1 loop dominates 

We first consider a typical scenario in which the dominant contributions come from the 
one- loop chargino and neutralino diagrams [32]. We set 

1^1 = IM2I = 2|Mi| = Mi^ = = Ms . (29) 

The common scale Ms and tan (3 are varied. For CP phases, we first consider the two 
values for $^ = 0° or 180° while taking vanishing CP phases for the gaugino mass and A 
parameters: $1.2,3 = $a = 0°. The remaining relevant parameters are fixed as 

Mh± = 300 GeV, 

Mq, = = = = = 0.5 TeV, 

\At,b,r^.\ = l TeV, |M3| = lTeV. (30) 

In Fig. [3l we show the one-loop contributions from the chargino and neutralino di- 
agrams (upper) and the two-loop contributions from the Barr-Zee graphs (lower) to the 
SUSY muon MDM as functions of Ms for several values of tan (3. The left frames are for 
$^ = 0° and the right ones for $^ = 180°. In both cases, the two- loop and one- loop contri- 
butions have same signs for most regions of Ms- The one-loop and two-loop contributions 
drops rapidly as Ms increases. The case of $^ = 0° gives the correct sign and, for example, 
we have (a^) susy x 10^° > 10 at Ms = 500 GeV when tan (3 > 20. For $^ = 180°, (a^) susy 
is negative and the considered regions of Ms are not allowed because data prefers a positive 
a^. Note that for tan /3 > 30, the tan /3-enhanced threshold corrections can turn the 6-quark 

9 




-0.5 



200 400 600 

Mg [ GeV ] 
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Figure 3: The one- (upper) and two-loop (lower) contributions to (a^)susY as functions 
of Ms ^29) for tan (3 = 2.5 , 5 , 10 , 20 , 30 , 40 , 50. In all frames, the smaller value of tan (3 
gives smaller \a^\. The other parameters are chosen as in Eq. ^3^). 



Yukawa coupling non-perturbative. This happens when Ms (or |/i|) is sufficiently large, as 
for the case at hand. It is shown by the termination of the curves in the right panels of 
Fig. [3l Our results are in good agreement with existing ones in literature. 

Figure H] shows the SUSY muon MDM (upper) and EDM (lower) as functions of 
$2 taking $i = 0° (left) and 180° (right). We have taken $^ = 0°, M5 = 250 GeV, 
tan (3 = 30, and the other parameters the same as in Fig. [31 In both the MDM and EDM, 
we observe that the dominant contribution is coming from the one-loop chargino diagrams. 
The subleading contribution from the neutralino diagrams is about 5 to 10 times smaller 
and likely has an opposite sign with respect to the dominant chargino contribution. The 
contributions from the Higgs-mediated Barr-Zee diagrams are negligible. Numerically, we 
have -0.8 < (a^)fusY x 10^° < 0.6 and \ {df,)^\ x 10^^ < 7ecm. We clearly see the (shifted) 
cosine and sine functional forms of the MDM and EDM, respectively, as to be anticipated 
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Figure 4: The MDM (upper) and EDM (lower) of the muon as functions of $2 taking 
$1 = 0° (left) and 180° (right) when = 0°, Ms = 250 GeV and tan/? = 30. The 
thick line is for the total MDM/EDM and the thin solid, dashed, dash-dotted lines are for 
the constituent contributions from the one-loop chargino, the one-loop neutralino, and the 
two-loop Barr-Zee diagrams, respectively. 



from the relations given in Eq. (ITT]) . In the upper frames, the horizontal band is the 
experimental l-a-allowed region, (30.7 ± 8.2) x lO^''; see Eq. (|3]). In the lower frames the 
region is overlayed with thick dots along the thick solid line. We note the chosen parameter 
set is compatible with the experimental data only for the non-trivial values of $2 around 
60° and 300°, resulting in large EDM of about ±5 x 10~^^ e cm, which can be easily observed 
once the projected sensitivity of 10~^^ can be achieved. 
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Figure 5: The muon SUSY MDM (left) and EDM (right) as functions of Mh, in the CPX 
scenario with tan/3 = 10, see Eq. [31\} . The lines are the same as in Fig. In each frame, 
the region left to the vertical line is excluded by data on T{1S) decay [50]. 



4.2 CPX scenario 

Next we consider the CPX scenario [48]: 



Mf 



/i| = 4MsusY 



M 



Da 



M 



Mi 



La ~ ^"Ea ' 
2 MsusY 



IM3I = 1 TeV. 



(31) 



Taking 
with $12 = 



A^, we fixed $a = ^3 



90°, and Mf 



M 



E2 



= 90°, MsusY = 0.5 TeV, IM2I = 2|Mi| = 100 GeV 
^susY- For our analysis, the most relevant feature 



of the scenario is that the combined searches of the four LEP collaborations reported two 
allowed regions where the lightest Higgs boson Hi can be very light for moderate values of 
3 < tan/3 < 10 [49]: 



Mh, < 10 GeV for 3 < tan/3 < 10, 

30 GeV < Mh, < 50 GeV for 3 < tan/3 < 10 . 



(32) 



On the other hand, a lower limit on the lightest Higgs boson, Mhi ^ 8 GeV, is available 
from the bottomonium decay T(IS') jHi [50]. Figure O shows (a^^) busy and dfj_ in the 
CPX scenario as functions of Mh^ taking tan/3 = 10. When Mh^ ^ 50 GeV, the one- 
loop contributions to (a^) susy are negligible compared to the Higgs-mediated two-loop 
contributions. The sign of (a/i)fusY pl^s (+) since it is dominated by the bottom-quark 
and tau-lepton loops mediated by Hi which is almost the CP odd state; see Eq. fl2S]) . 
However, it is still difficult to achieve (a^) susy x 10^° > + 10 only with a mostly CP-odd 
Higgs boson as light as ~ 8 GeV. For the EDM, the one- and two-loop contributions are 
comparable and tend to cancel each other. 
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4.3 An extreme scenario 

Finally, we consider a scenario in which the one-loop neutralino and chargino contribu- 
tions are suppressed while the two-loop Barr-Zee contributions dominate. This scenario is 
characterized by large tan (3, a light charged Higgs boson, very heavy smuons and muon 
sneutrinos, and very large and ,-1 parameters. Explicitly, we have chosen 



tan/3 = 50, Mh± = 0.2 TeV, 
= M^, = 10 TeV, 







- - - 1 TeV, 




\At\ = l TeV, 


IM2I 


= 2 |Mi| = 0.1 TeV, IM3I = 1 TeV, 






= 0°, 


$1,2 = 0° , 


(33) 



while varying 

^A,^ , $3 ; 1< |/"|/TeV < 12 , 1< I A^I/TeV < 50 . (34) 

where {Abrl = = \Ar\ and we have taken A^ = A-r. Note \At\ is fixed in this scenario 
and the results are almost independent of $1,2- 

Figure [6] shows the regions where (a^)susY = 30.2 ± 8.8 (blue) and (a^)susY = 
30.2±17.6 (blue+red) in the \AbT-\ and |/i| plane for several values of $3, taking ^a^^ = 180°. 
The unshaded regions are not theoretically allowed and we have M//^ < 100 GeV in the 
over-shaded regions (yellow), for example, in the upper-left corner in the upper- left frame. 
We found that Hi is always lighter than 100 GeV when $3 = 180° (lower-left) because the 
resummed threshold corrections modify the bottom-quark Yukawa coupling significantly 
in this case. We note the region with larger \Ahr\ is more preferred. Figure [7] shows the 
dependence on ^A^r of (a^)susY5 d^j,, and the masses of the lighter stau and the lightest 
Higgs boson. In the upper-left frame, the horizontal band is the experimental 1-a region 
of Aaj^^P. We found 150° < < 210° and ^a,, ~ 140° , 220° can make (a^)susY 

consistent with the experimental value for = 10 TeV and 10.5 TeV, respectively. In 
the upper-right frame, the 1-a region is overlayed with blank boxes along the dashed 
(|/i| = 10 TeV) and dash-dotted (|/i| = 10.5 TeV) lines. We have |ci^| x 10^^ < 120 e cm 
and dfj, x 10^^ ~ ±250 e cm for = 10 TeV and 10.5 TeV, respectively. We observe that 
the larger results in the lighter staus as shown in the lower- left frame. This leads to a 
larger (a^) busy as the dominant contribution in this case comes from the Higgs-mediated 
stau Barr-Zee graphs, as will be shown later. When = 5 TeV, Hi becomes lighter than 
100 GeV for ^a^^ ^ 80° and > 280°, as shown in the lower-right frame. When \fi\ > 10 
TeV, Mh, > 108 GeV. 

Figure [8] shows the dependence of {a^)svsY on (upper-left) and \Ahr\ (upper-right) 
for several values of \Abr\ and respectively, taking ^A^r = 180° and $3 = 0°. In the 
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Figure 6: Allowed region at the 1-a (blue) and 2-a (blue+red) level in the \Ai)r\ and \fi\ 
plane, see Eq. for ^3 = 0° (upper-left), 90° (upper-right), 180° (lower-left), and 270° 
(lower-right). We have taken ^Ai,^ = 180°. The unshaded regions are not theoretically 
allowed and M^^ < 100 GeV in the over-shaded regions (yellow). See Eq. for other 
parameters chosen. 



lower frames we also show the dependence of the mass of the lighter stau fi. Again we 
observe that large and \Abr\ can easily make (a^) susy consistent with the current Aa^^^ . 
Figure [9] shows various constituent as well as the total two-loop Barr-Zee contributions to 
the muon MDM (a^)'susY' functions of taking \AhT-\ = 40 TeV, ^A^r = 180° 
$3 = 0°. The thick line is for the total and the thin lines are for the constituent eight 
contributions (see Fig. [2]). For smaller the dominant contribution comes from the 
charged Higgs boson loop. This is because the Hi-H^-H~ couplings have loop-induced 
enhancement from large \Ai,\ and [51]. The possibility of such a significant contribution 
of the (photon-)Barr-Zee diagram with a (closed) charged Higgs boson loop has apparently 
not been noticed before. In fact, the particular diagram is typically not included in analyses 
at more or less the same level of numerical precision to ours presented here. As grows. 
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of the lighter stau (lower-left) and Hi (lower-right) for three values of fi: 5 TeV (solid), 10 



TeV (dashed), and 10.5 TeV ( dashed- dotted) . We have taken $3 
See Eq. [33) for other parameters chosen. 



0° and I A 



45 TeV. 



however, the contribution from the stau loops is enhanced and thus becomes dominant. 
The contribution from the sbottom loop is not significant because the resummed threshold 
corrections suppress the bottom-quark Yukawa coupling for the chosen parameter set. The 
scenario likely seems too contrived to some readers. However, we are presenting it here 
mainly to illustrate the significant roles of the various two-loop Barr-Zee contributions in 
some region of the parameter space. While the dominance of the latter group reduces as 
one moves away from the extreme corner of the parameter space, its significance maintains 
over a substantial region. For instance, we show the dependence of (a^)susY on M^^^^^ 
and \fi\ in Fig. [TOl i.e. the effect of bringing back the smuon mass from the heavy limit. In 
particular, the Higgs-mediated two-loop Barr-Zee contributions are shown to be actually 
dominate over the one-loop contributions even when the smuon mass parameters get down 
to as low as 1 TeV (left), for the full range of values (right). 
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Figure 8: Dependence o/(a^)susY on \fi\ (upper-left) and \ Ai,t-\ (upper-right) for several 
values of \Ai,t-\ and \fi\, respectively, taking ^a,,^ = 180° and $3 = 0°. The regions of 
\fi\ and are constrained by applying the experimental limit on the lightest stau mass 
as shown below. The lines are the same as in the upper frames. See Eq. ^33) for other 
parameters chosen. 



Before closing this section, we comment on the relation between the muon and electron 
EDMs. The most important one- loop contribution {di)^ and that from the two- loop {di}^ 
have somewhat different features. The Barr-Zee diagrams for the muon and the electron 
are identical except for the muon and the electron lines themselves. Hence, we have the 
robust relation 

(4)''= K/m^)(rf^)^. (35) 

For the case of (c?/)^, however, it is sensitive to the flavor dependence of the soft SUSY 
breaking terms. In most of the models on the origin of the soft SUSY breaking terms 
available in the literature, those of the first two generations are more or less the same. 
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TeV ] 

Figure 9: The Barr-Zee contributions to (a^)susY; ('2^)'susy? functions of |/i| taking 
\Ai)r\ = 40 TeV, ^A,,^ = 180° and $3 = 0°. Also shown are the constituent contribution 
from the eight types of diagrams. See Eq. [33^) for other parameters chosen. 



Explicitly, 
That does give us a relation 



{d,Y ~ {nie/ni,) {d,Y , (36) 



which is particularly sensitive to the CP phases of the Ae and A^ parameters. To summa- 
rize, one may consider special cases with [i] similar or universal soft terms (for electron 
and muon) or {ii) electron EDM dominated by the two-loop Barr-Zee graphs, de ~ (de)^ , 
possibly due to much heavier selectron masses. Depending on situations, d^ ~ {d^)^ + {dfj)^ 
or {d^)^ is very strongly constrained by the Thallium EDM as 

(i) |c/^| < 3 X lO^^^ecm ; (ii) |(c/^)^| < 3 x 10~^^ e cm , (37) 

where we used me/rrifj, = 1/213, c^ti — —585 ■ de and \dTi\ < 9 x 10~^^ecm [38]. Note that 
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dash-dotted (blue) lines are for the constituent contributions from the one-loop chargino, 
one-loop neutralino, and Higgs-mediated two-loop Barr-Zee diagrams. See Eq. for 
other parameters chosen. 



these limits are model dependent and the muon EDM should be measured independently 
of the Thallium EDM. 



5 Conclusions 



We have studied in detail various supersymmetric contributions to muon MDM and EDM, 
including one-loop chargino and neutralino diagrams and dominant two-loop Barr-Zee di- 
agrams. In general, the one-loop contributions dominate over the two-loop contributions, 
however there are interesting regions of the model parameter space where the two-loop 
Barr-Zee diagrams are the major source of contributions to muon MDM and/or EDM. 
The model parameter space is huge. It is not feasible for us to present and discuss here 
the numerical results for more than a few cases of interest. We try to pick cases that can 
illustrate the essential and interesting features, and leave it mostly to the readers to project 
onto the parameter space regions in between, or considered otherwise to be of special in- 
terest. We illustrate numerically 3 scenarios under various choices of soft parameters: (i) 
the one when the one-loop contributions dominate, (ii) the CPX in which the Barr-Zee 
dominates but the overall sizes of MDM and EDM are small, and (iii) a more exotic one in 
which the Barr-Zee dominates and the overall sizes of MDM and EDM are large. We have 
also shown interesting relations between the MDM and EDM. For the case when one-loop 
contributions dominate the MDM and EDM are just, respectively, the (shifted) cosines and 
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sines of the phase of the parameters involved. Existing experimentally preferred range of 
MDM already predicts an interesting range of EDM, which can be further tested in the 
future muon EDM experiments. For the case the MDM is dominated by two-loop Barr-Zee 
contributions the MDM and EDM can be connected by the relation in Eq. flMl) — a result 
of the more complicated Higgs sector phase structure. 

The CPX scenario may still allow a light Higgs boson after taking into account all the 
existing search limits. Potentially, the light Higgs boson could give large enhancement to 
muon MDM. However, after imposing the lower limit on Mh^ the resulting MDM is always 
less than 5 x 10"^", which is smaller than the experimentally favored value. 

The last scenario associated with large tan/3, large heavy smuons and muon 
sneutrinos, and large 1^6^-1 but light charged Higgs boson and stau is rather interesting. It 
suppresses the one-loop contributions but the two-loop contributions are large enough to 
explain the muon MDM data, which are dominated by charged Higgs boson and stau at 
smaller and larger respectively. For large enough the MDM data can be accom- 
modated easily. The particular interesting role of the charged Higgs boson escaped earlier 
studies. We have also illustrated that major features of the scenario persist over a region of 
the parameter space with milder conditions — in particular, more 'regular' smuon masses. 

In addition, we offer the following comments. 

1. We have included the threshold corrections to Yukawa couplings. In particular, the 
bottom-quark Yukawa can receive large corrections at large tan/5 with large IM3I and 

For tan/5 > 30 the threshold corrections can make the 6-quark Yukawa coupling 
turn non-perturbative when Ms (or is sufficiently large. 

2. The one-loop contributions to MDM and EDM can vary as shifted cosines and sines 
of the phase of the parameters. 

3. As shown in Fig. |U the prediction for EDM is of order 500 x 10"^^ e cm within the 
allowed range of MDM. It is about 2 — 3 orders of magnitude below the current limit, 
but will be within reach of future muon EDM experiments [41]. 

4. The CPX scenario may still allow a Hi as light as a few to tens of GeVs. It could be 
searched in the subsequent decay of the H2 — ^ HiHi, where H2 is the SM-like Higgs 
boson. The contribution of Hi to the muon EDM has a right sign but it may not be 
large enough to accommodate AaJ^^^ after taking account of the constraint from the 
bottomonium decay T(15') iHi. 

5. The last scenario that we studied is characterized by a very light stau, a light bino 
and wino, a light Mh-^, and a light charged Higgs boson. The predicted EDM is from 
to 100 X 10~^^ ecm. Experimental searches for this scenario at the LHC will be a 
lot of tau leptons in the final state because of lightness of stau. 
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6. The parameter space compatible with the AaJ^-^^ value required is generally extended 
by allowing CP phases. For example, in Fig. HI the CP-conserving cases ($2 = 0° 
and 180°) are excluded. 
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A CPsuperH Interface 

• Output: For output, part of auxiliary array RAUX_H is used. 

— The muon EDM in units of cm: 

RAUX_H(360) = dje = {djef^ + {d^/ef + {djef + {djef , (A.l) 

where the sub-contributions are 

RAUX_H(361) = {d^/eY^ , RAUX_H(362) = (rf^/e)*° 
RAUX_H(363) = (rf^/e)^ , RAUX_H(364) = (d^/e)^ . 

- The muon MDM: 

RAUX_H(380) = (a^) susy = (a;.)lusY + («m)Iusy + (aM)lusY + («M)fusY , (A.2) 
where the sub-contributions are 

RAUX_H(381) = (a^)lusY, RAUX_H(382) = (a^)|usY 
RAUX_H(383) = (a^)|usY , RAUX_H(384) = (a^)fusY ■ 

• IFLAG_H(19)=1 is used to print out EDM/MDM of the muon. Using run shell-script 
file distributed, the sample out obtained is 



The Electric EDM of muon in cm: 



d~E_mu/e [Total] : . 1288E-23 
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d"Ejnu/e[C,N,Gl,H] : O.OOOOE+00 -.6035E-24 O.OOOOE+00 0.1891E-23 



The SUSY MDM of muon: 



a_mu [Total] : . 1805E-09 
ajiu[C,N,Gl,H] : 0.1396E-09 0.2871E-10 O.OOOOE+00 0.1218E-10 
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